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We show that it is possible to associate univocally with each given solution of the time-dependent 
Schrodinger equation a particular phase flow ("quantum flow") of a non-autonomous dynamical 
system. This fact allows us to introduce a definition of chaos in quantum dynamics (quantum 

C ' chaos), which is based on the classical theory of chaos in dynamical systems. In such a way we can 

introduce quantities which may be appelled quantum Lyapunov exponents. Our approach applies to 
a non-relativistic quantum-mechanical system of n charged particles; in the present work numerical 
calculations are performed only for the hydrogen atom. In the computation of the trajectories we 
first neglect the spin contribution to chaos, then we consider the spin effects in quantum chaos. We 
show how the quantum Lyapunov exponents can be evaluated and give several numerical results 
which describe some properties found in the present approach. Although the system is very simple 

\£i ' and the classical counterpart is regular, the most non-stationary solutions of the corresponding 

t-H , Schrodinger equation are chaotic according to our definition. 
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> \ INTRODUCTION 
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y—{ • Around the turn of the century, Poincarc pointed out a property which is typical of most mechanical systems. 
| Although the evolution is deterministic, it shows an extreme sensitivity to initial conditions. Errors in these latter 
can grow up at an exponential rate, and thus in pratice, long-time predictions become impossible and the evolution 
appears irregular and random. In the modern literature this phenonenon is called "deterministic chaos" or also 
"classical chaos". [24| In the context of classical mechanics, it upsets the laplacian ideal of exact predictability, but, in 
the context of quantum mechanics, it rises even more profound problems. Indeed, in virtue of Bohr's correspondence 
principle, one expects that quantum mechanics must include or imitate, in some way, classical chaos. The direct 
y | analogue of this in quantum mechanics should be a sort of exponential sensitivity to the initial conditions in the 
solutions of the time-dependent Schrodinger equation. The difficulty to put in evidence such a phenomenon is well 
known [H, [2] and the reason is obvious. In quantum mechanics, instead of individual trajectories, we have a field: the 
(«j ■ wave- function yj(r,t), whose time-evolution is given by 



o 
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iP(r,t)=U(t,t )iP(r,t ), (1) 



where U(t 7 to) is a unitary operator [3[ on L 2 . Thus, according to L 2 -norm, two nearby initial wave-functions remain 
close as time evolves. Owing to this fact, which seems to rule out any possibility of a chaotic behaviour in quantum 
mechanics formally analogous to the classical one, the investigations on the subject have moved toward a different 
direction, that is in search of traces of classical chaos which persist in the quantum world. Although this research has 
obtained (see Ref. [3, [2] and references therein) interesting and far-reaching results, it looks less fundamental than the 
attempt of finding out a property of the quantum-mechanical time-evolution which is analogous to classical chaos. In 
the present paper, we show that such a property really exists. To put in evidence it, we give up, obviously, the too 
"rough" i 2 -norm to measure wave- function modifications and go to examine, so to say, "submicroscopical" details in 
the time-evolution (Q]). In short, these details are the "trajectories" r — r(t) determined by the condition of being 
tangent at any point r to the "velocity field" V(r, t) — J(r, t)/p(r, t) at time t and at point r, where J{r, t) and p(r, t) 
are, respectively, the probability density current and the probability density relative to the wave-function ip(r, t). We 
shall see that the set of these trajectories, distributed with the probability density p(r,t), can be considered as a 
particular phase flow (quantum flow) of an usual non-autonomous dynamical system, univocally determined by the 
wave-function ip(r,t) in (Q]). This latter fact permits us to introduce a definition of chaos for the time-evolution 
founded only on the basic concepts of the theory of chaos in dynamical systems. We call such a definition of chaos 
quantum chaos. As we shall see, our characterization of chaos in quantum dynamics implies also a sort of exponential 
sensitivity to the initial conditions in the Schrodinger equation and thus, it is the direct analogue of classical chaos. 
Notice that the emergence of quantum chaos, by "unraveling" the quantal evolution |T|) into "individual trajectories" is, 
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in a way, similar to the emergence of classical chaos in Halmitonian systems, by "unraveling" the time-evolution of the 
phase space probability distribution p(q,p, t) into individual particle trajectories. Indeed, since the equation (Liouville 
equation) satisfied by p(q,p,t) is linear, slight modifications in the initial phase space probability distribution never 
lead [|J to exponential divergences of p(q,p,t), although all or part of its "underlying" trajectories may be chaotic. 

It should be observed that our approach is formally similar to the one recently proposed [1, 0, [1] within the 
framework of the de Broglie-Bohm causal interpretation of quantum mechanics (Bohmian mechanics), however it 
presents the following novelties: 

- Our approach seems to be fuller. 

- We do not make recourse to any new interpretation of quantum mechanics. 

- We introduce suitable parameters (quantum Lyapunov exponents) to quantify quantum chaos. 

- We show a calculation of these parameters in a system simple but realistic (hydrogen atom). 

- We show intriguing relationships between quantum-mechanical quantities along a given trajectory and the 
Lyapunov exponents of this latter. 

Although all the numerical results presented (Sec. 5) concern the hydrogen atom, our characterization of quantum 
chaos and the relative implications (Sees. 2-4) are discussed in the general setting of a quantum system of n charged 
particles. 



THE QUANTUM FLOW 

In this section we show how it is possible to associate univocally with a given solution ip(r,t) of the initial value 
problem (IVP) for the Schrodinger equation a suitable phase flow of an usual non-autonomous dynamical system. 

Let us consider a system of n particles with masses m\,...,m n and charges ei,...,e n . The corresponding 
Schrodinger equation writes 



y — a, + v et6 ' 
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where r = (n, . . . , r n ), £ 1R 3 , i = 1,2, ... ,n. Each wave-function 

il>{r,t) = \p{r,t)] 1/2 j S(r ' t),h 
satisfying Eq. @ generates a probability density current 

J{r,t) = (j 1 (r,t),...,j n (r,t)), 



where 



and hence a "velocity field" 



where 



j t (r,t) : =P±A [ViS(r,*)] 



i = 1,2, . . . ,n, 



V(r,t) = { Vl (r,t),...,v n (r,t)). 



Vi(r,t) := 



z. = 1,2,..., n. 



(2) 



(3) 



(4) 



(5) 



(6) 



The vector field ([5]) is not defined at any point (r,t) E M 3n x Mt, where Mt denotes the time axis. In general, we 
can say that it is defined on some set 



V = R in x R t /S, 
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where 

S = {(r,t)e M 3n x M t : either tp(r,t) = (nodes) 

or ip(r,t) is not differentiable}. 

We call the set V extended phase space. We also need to consider the set 

V = {r G M 3n : (r, t) G V, for some t} 

and the set 

V(t) = {re M 3n : (r,i) G V, for t fixed}. 

We call the sets V and V(t) phase space and phase space at t, respectively. Let us now introduce in the extended 
phase space V the "direction field" (V(r,t),l). It is easily seen @, Sj that the curves in V, which at any point (r,i) 
have tangent directed as (V(r, t), 1) in that point, can be described by the parametric equations r = r(t),t = t, where 
r(t) satisfies the equation 

r(t)=V(r,t). (7) 

We name such curves integral curves. Notice that Eq. ([7|) defines what in the modern literature on chaos theory is 
referred to as non- autonomous dynamical system. For all we have said above, this system is determined univocally 
by the solution ip(r,t) of a given IVP for the Schrodinger equation 

iHd t ip(r,t) = Htp(r,t) 

^(r,0)=:Vto(r), Mr) € V{B), Uo{r)\\ = 1, [> 

where H is the Hamiltonian operator in (J2j) and T>(H) its domain. In association with the IVP (JH) we may also 
consider the following IVP 

r(t) = V(r,t) 

r(0)=:r o , r G V(0). [ J> 

The solution of ([9]), when represented in the extended phase space V, coincides with the integral curve which goes 
through the point (t"o,0). However, we prefer to think of the solution of (J9j) as the path (orbit) in the phase space 
V followed, as time develops, by a point starting from ro at t — 0. According to this representation, the solution 
of ^ will be called trajectory. Thus, the phase space V is the region (which can also coincide with the whole M 3n ) 
accessible to our trajectories and the phase space at t, V(t), is the subregion of V crossed by our trajectories at time 
t. [25| It should be observed that both the trajectory r = r(t) and the corresponding orbit are tangent at any point 
r to the velocity field V(r,t) at time t and at point r and that the orbit is the projection onto V of the integral 
curve. Obviously, in the spirit of the usual interpretation of quantum mechanics, the trajectories of the IVP ((9]) do not 
correspond to an effective motion of quantum particles. However, this fact does not have any relevance at all in our 
characterization of chaos in quantum dynamics. Indeed, as we shall see, what is most important to this aim is that the 
trajectories of the IVP ([9]) are formally analogous to trajectories of real point particles. Another consideration about 
our trajectories is that they must exist uniquely and globally in time, i.e. for all t G IRt- This property is, in fact, of 
basic importance to apply the theory of chaos in dynamical systems. Thus, we turn now our attention to it. First of 
all, we note that a prerequisite of this property is that the velocity field (|5|) exists uniquely and globally in time, i.e. 
the IVP © admits an unique and global solution. In order that it be true, all is required is [3j the self-adjointness 
of H, property which in our case is valid Consider now the main question: does the IVP © determined by ([8]) 
admit an unique and global solution? This question turns out to be more difficult, because the solution of ((9J may 
run into catastrophic events, such as reaching in finite time a point of the set S (the singularities of the velocity 
field). The question, which is also of basic importance in Bohmian mechanics, has been investigated recently [9j. The 
findings of the paper 0] are that, for a large class of potentials (including the Coulomb potential in ©) and under 
certain hypotheses (see for details), the IVP ([9]) admits an unique and global solution. This basic property will be 
assumed throughout the present paper. Now, in order to make more precise the definition of trajectory as well as the 
definition of quantum flow to be given below, we find it convenient to introduce the following mappings: 

a) $(-,ti, ri) : Mt — > V. This mapping represents the trajectory of the IVP ©, which goes through the point r\ 
at time t\. Thus Q(t, tx,rx) — r(t) represents the coordinate r at time t of the point pursuing this trajectory. 
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b) <!>(•, ti,-) : Mt x V(ti) — ► V. In the theory of dynamical systems, this mapping is called phase flow. In our 
picture, it represents the set of all trajectories of the IVP©, passing through V{t\) at time t%. 

c) $(£, ti, •) : V{ti) — > 'P(i). We call this mapping phase flow at t. 

The existence of these mappings is due B, [Toj to the assumption that the IVP (O admits an unique and global 
solution. This assumption implies also 0, the following properties: 

d) For t\,ti € Mt fixed, the mapping $(i2,ti,0 : T'(ti) — * ^(^2) is one-to-one, differentiable and with inverse 
differentiable (diffeormophism). 

c) •) = J (identity). 

f) $(i,t2 1 $(t2 1 ti 1 ri)) =$(i,ti, n) (chain rule) 

Now, let us introduce on V(0) a statistical ensemble £ of initial conditions (points), distributed according to 
the probability density p o (r(0)) = V ; o( T '(0))V'o( r (0)) = po(r(0)), where ipo( r ) is the wave-function appearing in the 
IVP ([5]), and let us consider the time-evolution of these points according to the IVP i.e. the phase flow $(•, •) 
defined on Mt x£. We call this special realization of the phase flow quantum flow (QF). Thus from one hand, the QF is 
formally similar to the phase flow of an ordinary dynamical system, on the other hand it has characteristic peculiarities. 
In fact, the velocity field V(r,t) in ([9]) depends on a specific solution of the IVP (JSj> and the initial conditions are 
distributed with a probability density, which is not uniform as in the phase flow of dynamical systems, but is given 
by po(r(0)). The fact that to each given solution ip(r, t) of © there corresponds an unique dynamical system ([7]) and 
hence an unique QF, gives us the opportunity of defining a chaotic behaviour in quantum dynamics. Indeed, we may 
say that the given ip(r,t) is chaotic when the corresponding QF presents trajectories which are chaotic, according to 
the usual definition of classical chaos. We call this characterization of chaos in quantum dynamics quantum chaos. 

As we shall show in Sec. 5, our definition of quantum chaos implies also a kind of extreme sensitivity to initial 
conditions, in the sense that a tiny modification in L 2 — norm in the initial condition ip(r,0) in (|8j gives rise to a 
dramatic modification in the QF pattern: each chaotic trajectory turns into another diverging exponentially from it. 
Strictly speaking, the definition of quantum chaos given as above is not fully rigorous and does not put in evidence 
the role of the ensemble £ . A rigorous and satisfactory definition of quantum chaos as well as a precise measure of it 
are postponed until Sec. 4, because first we need to describe some properties of the QF and also to recall the methods 
to quantify classical chaos. 



SOME PROPERTIES OF THE QUANTUM FLOW 



The properties of the QF are essentially those of an ordinary phase flow. However, due to peculiarity of the 
QF, some of them show themselves under new light and imply intriguing connections with quantum mechanical 
properties. Thus, besides mentioning the main properties of the QF, we dwell particularly upon these new aspects, 



giving also some proof, when this turns out to be necessary to make the paper clear and self-contained 27| . 



Property 1. Denoting by p(r(t),t) the probability density of points r(t) in the ensemble £ on V(t), the 
assumption, in the previous section, that the density of points in £ at t = on V(0) is p o (r(0)) = po(r(0)), 
implies HQ/ p(r(t),t) = ip*(r(i),i)ip(r(t),t) = p(r(t),t), Mt > and Vr(t) G V(t). 

Property 2. Let r$ be a given point ofV(0) and dro be the volume element around this point r 2&]. Let r(t) and 
dr(t) be the corresponding quantity tranformed by the QF at t > 0, i.e. r(t) = $(t, ro), dr(t) — &(t,dr ). Then, the 
probability measure p{r(t),f) dr(t) is invariant under the QF, i.e. 

Po(r Q )dr = p(r(t),t)dr(t), Vi > 0. (10) 



Proof. Think of the QF at t (cf. points c) and d) in Sec. 2) as the coordinate transformation 

S(t,r(0)) = r(t) (11) 

from V(0) to V(t), t > 0. We have 

dr(t) = J(t, ro) dr , (12) 
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where J(t,ro) denotes the Jacobian at ro of the transformation (fTTj) . This Jacobian satisfies [ll[ the equation 

j= (j2 Vi-vA J, (13) 

where the quantities Vi are given by ©. Integrate (H2J) with the initial condition J(0, ro) = /, and integrate the 
continuity equation for p in the Lagrangian form, i.e., 

^+p(Ev,-^=0, (14) 

with the initial condition po(?"o)- We find 

p{r(t),t)J(t,r ) = Po {r ), Vi > 0, 
which, in virtue of (fL2|) . yields the relationship (fT0|) .D 



Property 3. The growth rate under the QF of the size of the volume element dr(t) around the point r(t) is given 

by 

,n!rU) =dr(t) [J2 Vi-Vi(r(t),t) ) , Vt>0. (15) 



df 



The proof is an easy consequence of Property 2 and of {33]). The integration of Eq. (fT5")) over a finite volume 
V(t) C -P(t) yields 



eft 



/ V ( •«<(»•(<),«) j dr(i). 



Since (see Eqs. ([51315]) ) 



n n 1 

£v i -t> i (r(t),t)=X; — (A,S(r(f),i)), 



, mi 



we find that in a given region f2 C T 5 , the QF is conservative or nonconsevative, according to whether the value of 
S"=i W~ (AjS'(r(t), t)) in O is equal to zero or different from zero, respectively. 

We have to recall now the basic elements to define and quantify chaos in dynamical systems, i.e. the definition and 
the main properties of the Lyapunov exponents (LEs) of a given trajectory. We do that going in some details, because 
these are necessary to put in evidence the interesting connections mentioned at the beginning of this section as well 
as to explain part of the numerical results shown in Sec. 5. 



According to Refs. [12|, [13J, the LEs of the trajectory <!>(•, r ) starting from r at t = and belonging to our QF, 



arc defined by the following limit 



A(r ,£ ):= lim - In \\U ro (t) £ || , (16) 

t— >+oo t 



as function of £ , where || • || is the Euclidean norm, £ is an arbitrary starting vector of M 3n , ||£ |j = 1, and U ro (t) 
is an in x 3n invertible matrix (flow matrix) . This matrix is the solution of the IVP 



•(■.ro)Ju = r™ (17) 



U ro (0)=I, 



obtained by linearization of ([9]). Rigorous results on the existence of the limit (| 16[) and on the properties of the 
corresponding LEs can be found in Refs. 12, [l3| and references therein. It is not easy to show (it would be also 



beyond the scope of the present paper) that the findings in [12j, |13|, proved under suitable hypotheses, can be taken 
over to the dynamical system ([?])• However, on the basis of a strong numerical evidence found by us in the case 
3n = 3, it seems reasonable to extend to ([7]) the validity of what follows: 
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i) Depending on the choice of £ , the limit (| 16[) takes on 3n (not necessarily distinct) values: 

Ai(r ) > A 2 (r ) > . . > A 3 „(r ). 

ii) If £ is chosen arbitrarily, the value of the limit (fT6|) is Ai(ro). 

iii) Consider the symmetric and positive matrix 

1/2* 



(t):=(u ro (t)U ro (t)) , (18) 

where U denotes the transpose of U. Denote by Hi{t, r ) > . . . > ^3n(^, r ) its eigenvalues. Let 

Xi{t,r ) := bxm(t, r ), i = l,...,3n. (19) 

Then, 

lim Ai(t,7- ) = Ai(r ). □ 

t — >+oc 

The meaning of the LEs is put well in evidence by the following geometric argument. Consider the singular value 
decomposition [3| of the matrix U rQ (t), i.e. 

U ro (t)L ro {t)=G ro (t)F ro (t), t>0, (20) 

where L and G are orthogonal matrices and F is a diagonal matrix with elements o i it, 7*o) > 02(^,7*0) > ... > 
C3n(i, ro) > 0. Take account of the fact that ([20]) implies 

so that the columns of L are the orthonormal eigenvectors ei(t, Tq), . . . , e 3n (t, Tq) of M^-q (t) and 

ai(t,r ) = exp[Xi(t,r )t], i = l,...,3ra. (21) 

Then, 

U ro (t)ei(t,ro) = exp[A i (t,r )t]e-(t,r ), i = 1, . ..,3n, (22) 

where e£(i,T*o) are the columns of G. Let now dro be an infinitesimal parallelepiped, with axes £iei(t,ro), i = 
1, . . . , in, which at t — is placed at rg. Since the Ei are infinitesimal, 

U ro {t)e l e l {t,r ) ~ ^>{t,r + e t e,(t,r )) - $(t,r ) , 

i = 1, . . . , 3n , 

and thus, in virtue of (|22| . the phase flow at t > transforms <iro into another infinitesimal parallelepiped dr(t), 
which presents a different orientation in lR 3n and the iih axis, i = 1, ...,3n, grown or shrunk according as the 
corresponding Ai(i,T*o) is positive or negative, respectively. The value of \i(t, ro) gives the rate at time t of the 
exponential growingor shrinking of this axis. The meaning of A 2 (ro), i = 1, . . . , 3n, is then obtained accordingly in 
the limit t — ► +00 [3fJ . We can now prove the following property: 

Property 4. Let dro and dr{t) be volume elements the same as in Property 2. Suppose dro is a parallepiped. 
Then, 



dr(t) = dr exp 



3n 



vt > o . (23) 



Proof. In virtue of (|22[) . each axis 

d£i = £\ei{t,r ) + . . . + el n e 3n (t,r ), i = l,...,3n, 
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of dro is transformed according to 

Cr (*)d£i = d£< = 4exp[Ai(t,r )t]e' 1 (t,r ) + ...+ 
4« exp[A 3 „(t, r )t]e3 n (t, r ). 

Taking account of this fact and recalling that the 3n-dimensional volume spanned by 3n linearly independent vectors 
arranged in an 3n x 3n matrix A is | det(A)|, the proof is completed by a straightforward calculation. □ 



The relationship (|23[) entails two more properties. 
Property 5. 

p(r(t),t) = p {r ) exp 



3 /j 



Vt>0, (24) 



where the points r(t) and r$ are the same as in Property 2 and p(r(i), t) = ip*(r(t),t)ip(r(t),t). 
Property 6. 

n _^ 3n 

V — (AiS(r(t),t)) = Y][\i(t,r ) + \i(t,r )t\ , 

1=1 1 i=l 

Vt > , (25) 
where £ is the phase ofip(r, t) in the polar form and the points Vq, r(t) are the same as in Property 2. 

The proof of (f2"4"|) is a straightforward consequence of (f2"3"| and of Property 2. To prove (|23|) . it is enough to 
differentiate with respect to t Eq. (f23|) and then to compare the result with p5|) and to recall ([4]) and ([6]). These two 
latter properties seems to us somewhat intriguing, for the following reason. The quantities A^(t, t*o) are functions of 
t and hence of the point r(t) = $(t,ro). Therefore, they may be looked upon as "classical dynamical variables" of 
the point pursuing the the trajectory <!>(•, ro). In virtue of the relationships and (|25p. these variables are strictly 
correlated to the time-evolution of ip{r , t) along the given trajectory. 



QUANTUM LYAPUNOV EXPONENTS 

In this section, we give the precise definition of chaotic solution for the IVP ^ and at the same time we give the 
tool to measure the corresponding degree of quantum chaos. We get this tool by considering a properly weighted 
average of the LEs of all trajectories in the associated QF. 

Definition 1. Let ip(r,t) be the solution of the IVP 10). Recall that po(r(0)) is the corresponding probability 
density of points in our statistical ensemble E at t = 0. Then, the quantities 

Ai(0):=/ Ai(r(0))p (r(0))dr(0), i = 1, 2, . . . ,3n, (26) 
Jv{o) 

where Ai(r(0)) is the ith LE of the trajectory of (0) starting from r(0) at t = 0, are called "quantum Lyapunov 
exponents" (QLE) of the given ^(r,t). For the sake of simplicity hereafter we suppose the usual normalization 
N = J v{0) p o (r{0))dr{0) = 1; if N ^ 1, we must divide by N the r.h.s. of eq. (26). 

Property 7. The QLEs are invariant under the QF, i.e., 

A;(0) =Ai{t) := [ X l (r{t),t)p{r(t),t)dr(t) , 
Jv(t) 

i = l,2,...,3n, Vt > , 



where r(t) = <!>(t, r(0)) and Aj(r(t), t) is the LE of the trajectory of (0) starting from r(t) at t. 
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Proof. As in the proof of Property 2, think of the QF at t as the coordinate transformation Q(t, r(0)) = r(t), from 
V(0) to V(t). We find 

Ai(t) := / Ai((r(t),t))p((r(t),t))dr(t) 
= / Ai(*(i,r(0)),*)po(r(0))dr(0) 

J-P(O) 

= / Ai(r(0))po(r(0))dr(0) 
Jv(o) 

=: Ai(0), i = l,2,...,3n, Vt > 0, 

where we have employed (jTOJ) and the property Xi(^(t,r(0)),t) = Aj(r(0)) which follows from the very definition of 
LE.D 

According to the above property, the QLEs are parameters whose values are univocally correlated to peculiar 
features of a given quantal evolution, i.e. the chaotic behaviour of the trajectories in the associated QF. Therefore, 
in correspondence with the definition of classical chaos, we can introduce the following definition: 

Definition 2. A given solution of the IVP {HP is called "chaotic" if the leading QLE Ai [31] is positive. The value 
of Ai is assumed as a measure of the corresponding quantum chaos. 

Notice that employing only Ai in the definition above corresponds to the fact that the quantity of main interest, 
to quantify chaos in a dynamical system, is the leading LE. Another quantity is the so-called KS entropy 0, EEI 
which requires all the positive LEs. In principle, an analogous parameter employing all the positive QLEs, could be 
introduced also here, but we do not find it necessary to discuss this point further in the present paper. 



NUMERICAL RESULTS 



Since the leading QLE is a precise indicator of our characterization of quantum chaos, most of the present section 
is devoted to discuss in some detail, through a numerical example, how it is possible to calculate it. We also show 
numerical results to illustrate some properties of the QF and the complex variety of the underlying trajectories. As 
is evident from the definition (|26p . the computation of Ai is, in general, very difficult. Indeed, once assigned a given 
solution [il] ip(r,t) of the IVP flHJ, this calculation requires what follows: 

1) Numerical integration of the IVP associated with the given wave-packet ip(r,t). 

2) Calculation of the LE Ai(r(0)) for all trajectories in the QF relative to ip(r,t). 

3) Evaluation of the integral (|26)l . in the case i = 1. 

As we shall see, each point above presents peculiar numerical problems. The numerical example shown below is 
relative to the calculation of Ai , in the case of the wave-packet 

i>(r,t) = ^100 (r)e5* + ip 200 (r)e** +^211 (r)e**) 

= p 1 / 2 (r,i)e lS ( r ^\ (27) 

where ip n lm(f) denote eigenf unctions of the hydrogen atom [33j ]. The chosen wave- function is the most simple su- 
perposition of the hydrogen atom eigenstates entailing, as we are going to show, the presence of quantum chaos and 
non-vanishing QLEs. As concerns point 1), we employ an adaptive fourth-order Runge-Kutta method, evaluating nu- 
merically the velocity field v(r, t) = VS(r, t) appearing in ([9]). We find that a considerable portion of the trajectories 
are chaotic, i.e. Ai(r(0)) > 0. This fact implies the well known difficulty [lj| about the integration of chaotic systems: 
numerical errors propagate exponentially and the trajectories after a certain time, depending on the accuracy of the 
algorithm employed, become meaningless. This phenomenon is amplified if, as in our case (see below), there is no 
attractor. However, rather than in trajectories, we are interested in the corresponding LEs, and these latter (see 
Fig. I) are stable to within 10%. 
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As concerns point 2), we should evaluate only Ai(r(Q)), the largest of the three LEs which appear in our example. 
However, since we want to check the relationship (|24|) and also for reasons (see below) of numerical control, actually, 
we calculate all three LEs. It is well known [lg, ll7[ that, contrarily to the evaluation of the only leading LE, the 
calculation of all of them present numerical difficulties and special techniques are called for. In the present paper we 
have employed the algorithm introduced in Ref. [17| (BGGS algorithm) and described also in Ref. [161 ]. 

This algorithm is based on formula ([TB|), when one is interested only in the evaluation of the leading LE, otherwise 
it employs the more general formula 



Ai(ro) + • • • + A p (r ) 
1 



km j\nVoF{U ro (t)^... U ro (t)£ p }, 

— . -4- rv~i t 



l<p<3n, (28) 

where £ 1; . . . , £ p are orthonormalized starting vectors arbitrarily chosen, Vol p {-} denotes the volume of the paral- 
lelepiped defined by the vectors U ro (t)^ 1 , . . . , U ro (t)£ p and r<j is the specific point from which the trajectory starts at 
t = 0. The peculiarity of the BGGS algorithm is that of making use of successive orthonormalization procedures in 
order to avoid well-known 16|, IjJ problems of overflow and of "alignment" of the starting vectors. With this trick, 



the BGGS algorithm turns out to be very effective and preferable to the evaluation of the LEs directly as logarithms 
of the eigenvalues of the matrix (|18p. Notice that the BGGS algorithm yields, in pratice, approximations Aj(£, ro), 
1 = 1,..., 3n, as functions of t for the LEs. In comparison with the analogous quantities (fT§|. the BGGS approxima- 
tions are, as we said, more reliable numerically. On the other hand, contrarily to (|19p . they depend on the starting 
vectors £ l5 . . . ,£ p , and thus, strictly speaking, they cannot be considered as the "dynamical variables" which enter 
into the relationships (|23p - (l25|) . However, in spite of the fact that the quantities (fl9|) and the BGGS approximations 
coincide only in the limit t — > +00, the respective sums of all them coincide at each t. This latter property may be 
seen comparing Eq. ([23| with the analogous one obtained by truncating at t Eq. (|28|) in the case p = 3n. In virtue of 
this coincidence, we shall employ the quantities Xi(t,ro) obtained by the BGGS algorithm to check relationship (|24| . 
The numerical difficulty of point 2) is that the wave-packet l|27p . nothwithstanding its semplicity, gives rise to a 
complex scenario of trajectories. In fact, we find trajectories whose behaviour passes, through all the intermediate 
degrees, from the regular one to the chaotic one including also the so-called "intermittency" [2j. In this latter case, 
the trajectory switches back and forth between two qualitatively different behaviours, which in our case are the mildly 
chaotic one and the chaotic one. To put in evidence these different regimes, we show a typical regular trajectory 
(Fig. 2a), a typical "intermittent" chaotic trajectory (Fig. 2b) and a typical chaotic trajectory (Fig. 2c), and for each 
of them, together with the LE Ai(i, ro) [34| as function of t, we display the quantities 

\W __Li W u r{t k ){tk+\,tk)i{tk)\\ . . 

where k = 0,1,2,..., r(t k ) = $(i fe ,r ), t k = k At and £(*k+i) = U r ( tk ){tk+i-,t k )£,(t k ). The symbol U r{tk) (t k+ll t k ) 
denotes the flow matrix solution of the IVP Q17p with initial condition at t = tk, U r ( tk ){t k ,t k ) = /. In other words, 
the A< fe) measure the rate of the exponential divergence or convergence, in the time interval (tfe+i,tfe), of the two 
trajectories: $(•, r ) (the fiducial trajectory) and $(•, r +e^ ), where e is infinitesimal and £ is an arbitrary starting 
vector, ||£o|| = 1. Furthermore, recalling the chain rule for the flow matrix (cf. point f) in Sec. 2), the definition (fTB|) 
and point ii) in Sec. 3, we have 

1 1 N ^ 

M(t N ,r ) := — ln\\U ro (t N )U = „ £ A«. 

N k=0 

Strictly speaking, Ai(tjV)T"o) shows a dipendence on £ , which disappears in the limit t — > +00. Here, we do not find 
it necessary to investigate such a dependence, because in our case the LEs present a relative error of order 10%. The 
above mentioned complex scenario of trajectories has forced us to find out a numerical device (see below) capable of 
giving a reliable estimation of 

Ai(r(0))= lim Ai(t,r(0)), (30) 

' + 00 

over the whole range of the chaotic behaviour. The point 3) in our example reduces to the evaluation of 

Ax = / A!(r(0))p (r(0))dr(0), (31) 
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where p (r(0)) is the probability density in (|27|) at t = 0, r(0) varies on 1Z C 7^(0) and 7?. = {(x.y, z) : — 9 < x < 
9, — 6 < y < 6, — 8 < z < 8} is a region such that the density po^O)) is negligible outside. The numerical difficulty 
here is that the values of Ai(r(0)) are given, as we said, with a relative error of order 10%. For this reason, we have 
found it convenient to employ an adaptive Monte-Carlo method and to perform four integration runs with different 
accuracy parameters, 35[ to see how much the errors over Ai(r(0)) affect the value of Ai. 



We pass now to describe how the actual evaluation of (f3"Tj) is carried out. Chosen an initial condition r(0) € K by 
the Monte-Carlo routine, the integration of the IVP ([9]) and the BGGS algorithm start in simultaneously. After a 
time, say t\, large enough, we apply the numerical device previously mentioned to estimate the limit (|30|). Actually, 
this numerical device consists in two different convergence criteria. The first one is conceived to recognize regular 
trajectories and it is as follows. We divide the interval I\ = [0.4ti,ti] in five sub-intervals of equal lenght and 
compute a linear least squares approximation (LLSA) to Ai(i,r(0)) in each sub-interval. [36| If all the LLSAs are 



decreasing and the minimum of Ai(i,r(0)) on I\ is smaller than a fixed threshold, [37j the corresponding trajectory 
is considered regular and its LE Ai(r(0)) in (|30| is estimated at zero (Fig. 3). As concerns the second convergence 
criterion, conceived for chaotic trajectories, we compute a LLSA of the function Ai(t, r(0))t on 1%. If the root-mean- 
square-error in the LLSA is sufficiently small, the angular coefficient of the fitting straight line is assumed as value 
of Ai(r(0)) (Fig. 4). The integration of the IVP © and the BGGS algorithm stop, if one of the two convergence 
criteria is satisfied, otherwise they go on up to a certain time ti > t\ and the two criteria are again applied on the 
new interval [0.4^2, £2], and so on, until one of the two criteria is met. Then the procedure starts again with a new 
initial condition. 

In order to evaluate the QLE (j3"T|) with a relative accuracy of 10%, the Monte-Carlo routine has employed about 1700 
initial conditions and the corresponding CPU-time in a double-precision calculation on a workstation, has been about 
60 hours. As we said, we have performed four such runs, varying the accuracy parameters, and the correspondent 
values found for Ai are the following: 9.048 • 10~ 3 , 9.326 • 10" 3 , 8.465 • 10" 3 , 9.048 • 10" 3 . 

Nothwithstanding the scenario of the trajectories pertaining to the wave-packet (|27|) is complex, we have found 
that our integration domain 1Z may be subdivided into three sufficiently disjoint regions, according to the Ai(r(0)) 
values. These three regions (Fig. 5) are: the regular one (Ai(r(0)) = 0), the mildly chaotic one (0 < Ai(r(0)) < 10 -2 ) 
and the chaotic one (Ai(r(0)) > 10~ 2 ). On the basis of the values found for Ai and considering also that the most 
chaotic trajectories have, in our units, LEs of order 0.04, we can say that the wave-packet ([27)) presents a good 
degree of quantum chaos. Although we have investigated in detail only the wave-packet (|27p . we have the numerical 
evidence, according to what observed in Ref. (Hj for a different system, that there are chaotic trajectories in the QFs 
corresponding to most wave-packets relative to the hydrogen atom. Thus, even if the classical counterpart of the 
hydrogen atom is not chaotic, our characterization of quantum chaos seems to spread through the corresponding 
quantum dynamics. However, there are exceptions, such as the following: 

i) Spherically symmetric wave-packets ip{r,t). 

ii) Wave-packets of the kind ij}{ r , t) = Ek a fc/fc( r ) e lEkt where a k , fk € M, and the energies E k are commensurate 

iii) Stationary wave-packets ip(x,y, z)e~ lEt , where (p(x,y,z) — ip* (x, —y, z). For example, ip(r,6,(j>,t) — 
e lEt Ek <Xkfk(r,0) e im "<t>, where a k J k eM,m k =0,±1,.... 

The wave-packets i), ii) always give rise to regular QFs, whereas the wave-packet iii) only in some situations. Although 
the explanation of these facts is not very difficult, it cannot be given in few words, and it will be shown within a 
more general framework elsewhere. Here, we confine ourselves to note that iii) includes, in particular, the stationary 
solutions R n i(r)Yi m (9,(j))e~ lErit of the time-dependent Schrodinger equation for the hydrogen atom. In this latter 
case, all orbits are circles centred on the z-axis and parallel to the (x, ?/)-plane. Obviously the corresponding Ai is 
zero. 

Fig. 6 shows the check of the relationship (f2"4"| . Although this relationship holds, in principle, for any t > 0, in 
practice, its validity is resticted to a finite time interval < t < T, where T depends on the numerical accuracy. In 
our case T is very large (of order 10 6 a.u.). In Fig. 6, to put in evidence discrepancies in a short time interval, we 
have reduced the accuracy in the calculation of the LEs. 

Fig. 7 illustrates that 

3 

lim VAi(t,r(0)) =0, Vr(0) e TZ. (32) 

i=l 

In virtue of the relationship (HM)) . limit (13"2")) means that our trajectories do not reach the nodes of the wave-packet P?]) 
and that the density p(r(t),t) keeps bounded. Actually (see Fig. 6), the density p(r(t),t) along a given trajectory 
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shows bounded oscillations. Thus, recalling Property 2, we can also say that our QF, as concerns the variation of the 
size of volume elements, has not any definite character, but it is expected to be conservative on average. We have 
employed sistematically the limit relationship (|3"2"|) to check the accuracy of the LE calculation. 

Finally, Fig. 8 illustrates the fact that the pattern of a chaotic QF is very sensitive to tiny modifications in L 2 -norm 
of the initial conditions in the IVP © ■ To give a clear idea of this phenomenon, Fig. 8 displays only the components 
z(t) (as functions of t) of the points pursuing two typical chaotic trajectories, both starting at t = from the same 
initial point, but pertaining to two diffcrnt QFs. The first QF is relative to the IVP |J5]) with initial condition i/jo(r), 
obtained by setting t = in (f27|) . The second QF is relative to the IVP ([8]) with initial condition ip' a (r), obtained from 
tpo(r) by a tiny modification of the coefficients of the corresponding linear combination, such that ||^>o — tp'oW = 10~ 4 . 



SPIN CONTRIBUTIONS FOR THE QLE 



Whenener we deal with a spinning particle — as is the electron of the hydrogen atom — we must take into account 
also the presence of the intrinsic angular momentum. The presence of the spin in the Schrodinger theory does not 
affect the stationary energy levels and the wave-equation to be used is, as usual, Eq. ((2]). Nevertheless, we have 
to consider the spin, not only to obtain the allowed stationary states of atoms (because of the Pauli principle and 
the "exchange interaction"), but also to write the complete expression of the conserved probability current for the 
Schrodinger equation. The complete Schrodinger current can be got [3, 2, 2^, [2l|: a) taking the non-relativistic 



limit of the Dirac current; b) considering the Pauli current in the particular case of a spin-eigenstate wave- function; 
c) directly from the Schrodinger equation, by considering the privileged direction which the spin introduces in the 
space. Let us remember that any solution of the Schrodinger equation is always a spin eigenstate along some space 
direction, because the Hamiltonian —ft 2 A /2m + U does not contain any spin term. Then the "complete" form for 
Schrodinger wave- function for a spinning particle actually is \fp^ lS ^ n X where \ is a constant 2-components Pauli 



spinor. Following Landau 18[ we can write the conserved Schrodinger current as follows: 



i = -VS+-^-Vpx S . (33) 
m 2m 

where s is the constant unit versor in the direction of the spin: s = X o'X/x Xi 47 being the usual Pauli 2x2 vector 
matrix. [38| When the polarization direction is chosen parallel to the z— axis, i.e., s — (0; 0; s z ), s z = ±1, we have 
for the spin-part of the current: 

— Vp xs = s z (V yP ; -V x p; 0) . (34) 
2m 

The current in Eq. (|33[) appears as a sum of the translational, "newtonian", part pVS/m which, at the classical limit 
(ft — > 0), is parallel to the classical impulse (equal to ft times the gradient of the action); and of a "non- newtonian" , 
rotational part due to the spin, Vp x s/2m, which vanishes only at the classical limit, i.e. for spinless bodies, but 
not in the non-relativistic limit, i.e. for a small momentum p. Let us stress that for any solution of the Schrodinger 
equation the spin part of the current is conserved because the curl of a divergence identically vanishes: 

V • [Vp x s] = V ■ [V x (pa)] = ; 

as a consequence, if pVS/m is conserved, quantity pVS/m + ftV x (pa) /2m is conserved as well: 

p ft 



dtp + V • (—VS) =d t p + V ■( ^-VS + —Vp xs =0 

V m / \ m 2m I 



Let us also remark that the spin part of the conserved current results to be of the same order of the scalar-newtonian 
term] actually, in the case of real wave-functions, the spin part is the non- vanishing one in the probability current. (39j 
Let us finally recall that also in Bohm Mechanics the inclusion of the spin part in the non-relativistic velocity field is 
requested for spinning systems 0, 0|- For a n— particles system Eq. (|4]) generalizes as follows in the presence of spin: 

3i(r, t) := V,5(r, t) + JL V lP (r, t) x Si . (35) 

mi 2mi 

while the velocity field Jj(r, t)/p, which for a spinless particle consists of the only term ViS(r,t)/m^ now becomes 

Vi (r, t) := — V,5(r, t) + ^—-V t p(r, t) x Si . (36) 

mi 2mip(r,t) 
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As it is easy to check, all the properties of the QF proved in Sec. 3 still equally hold in the presence of spin. We have 
performed several numerical calculations of the chaotic quantities for the hydrogen atom considering the spin: that is, 
by integrating the above "complete" velocity field. We get results quite similar to the ones obtained in the absence of 
spin and discussed in the previous section, but the spin part of the quantum flow seems to increase slightly the chaos 
with respect to the spinless case. As an example, for the wave-packet 



with spin projection s z — +1, we have obtained (after four runs, varying the accuracy parameters, etc.) the following 
value for the QLE: Ai = 1.022 • 1(T 2 . 



In this paper, we have proposed a precise definition of chaos in quantum dynamics. This definition is based only 
on basic concepts and tools of the theory of chaos in dynamical systems and, contrarily to previous papers, it does 
require any new interpretation of quantum mechanics. The characterization of chaos in quantum dynamics introduced 
by us may be looked upon, in a sense, as the direct analogue of classical chaos and thus, it attempts to reduce the gap 
between classical chaos and quantum chaos according to the definition in most current literature. The present work 
stimulates further investigations mainly on two fundamental problems. The first one concerns the correlation with 
classical chaos. We expect that, if we adopt the criterion for the classical limit described in Ref. quantum chaos 
here introduced matches with classical chaos. The other problem is that of seeing if the values of the QLEs may be 
related to observable quantal quantities. We want to investigate these two problems in a next paper and work in this 
direction is already in progress. 

The authors would like to thank G.G.N. Angilella, M. Baldo, M. Baranger, S.P. Goldman, V. Latora, F. Raciti, A. 
Rapisarda and E. Recami for many helpful discussions. This work has been supported by INFN, Sezione di Catania, 
and by MURST. 
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FIG. 1: Histogram relative to the distribution of the values found for Ai(ro) over 100 trajectories, all starting from the same 
initial point r but evaluated with slightly different accuracy. N is the number of trajectories, (A) is the expectation value and 
a is the standard deviation. 
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FIG. 2: Figures (a)-(c) are relative, respectively, to a typical (regular, "intermittent" chaotic, chaotic) trajectory. The solid 
curve represents the LE Ai(£, ro) as function of t and the dots represent quantities (29). 
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FIG. 3: Convergence criterion for a typical regular trajectory starting from ro- The solid curve represents the LE \i(t, ro) as 
a function of t. The dashed straigth lines represent the LLSAs. 
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FIG. 4: Convergence criterion for a typical chaotic trajectory starting from ro. The solid curve represents \i(t,ro)t as a 
function of t. The dashed straight line represent the LLSA. 
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FIG. 5: Projection onto the (y, 2)-plane of all the initial conditions in the region {(x, y,z) STZ : —0.5 < x < 0.5} and relative 
to the trajectories evaluated in the four runs to calculate Ai. Three regions of different degree of chaos are put in evidence. 
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FIG. 6: Check of Eq. (|24[) . The solid curve represents the probability density p(r(t),t) as given by the wave-packet (|27[1 along 
a typical chaotic trajectory. The dashed curve represents p(r(t),t) along the same trajectory but evaluated through (|24l) . 
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FIG. 7: Check of the limit relationship (|32|l . The solid curves represent the three LEs: Xi(t, ro), M(t, ro), \3(t, ro) as functions 
of t, of a typical chaotic trajectory starting from ro. The dashed curve represents Y^r=i ^i(t> r o)- 
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FIG. 8: Dependence of the QF pattern on the initial condition in the IVP (JSj) . 



